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Abstract 

A long Josephson junction in a constant external magnetic field and in the 
presence of a dc bias current is investigated. It is shown that the system, simu- 
lated by the sine-Gorgon equation, "remembers" a rapidly damping initial per- 
turbation and final asymptotic states are determined exactly with this perturba- 
tion. Numerical solving of the boundary sine-Gordon problem and calculations 
of Lyapunov indices show that this system has a memory even when it is in a 
state of dynamical chaos, i.e., dynamical chaos does not destroy initial informa- 
tion having a character of rapidly damping perturbation. 
PACS number(s): 74.50+r, 05.45. +b 

Dynamical chaos is one of the most interesting phenomena in the 
theory of Josephson junctions. 1-10 This phenomenon is not only of the- 
oretical importance but also of practical importance, because many de- 
vices are founded on Josephson junctions, in particular, superconducting 
quantum interference devices (SQUID's). 11 Dynamical chaos in these de- 
vices is another source of noise. Furthemore, a long Josephson junction 
(LJJ) serves as a very good system for studying nonlinear phenomena 
such as an exitation of fluxons and antifluxons, their propagation, in- 
teraction, scattering, and breakup. Investigations of the last few years 
showed that a LJJ detects deeper characteristics that had seemed. Even 
in the simplest case, when a bias current and an external oscillating 
field are absent, the presence only of a constant external magnetic field 
leads to the most interesting phenomenon connected with the selection 
of the solution of the stationary Ferrell-Prange equation. The fact is 
that this equation has not only provided one solution by given boundary 
conditions; the number of these solutions increases by the strength of 
the external magnetic field and the total length of the junction. 11 
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Recently we have shown that the selection of a solution is carried out 
with the form of a small and rapidly damping initial perturbation in time 
in the nonstationary sine-Gordon equation and what is more surprising 
an asymptotic solution of this equation coincides with one of the stable 
solutions of the stationary Ferrell-Prange equation. Two circumstances 
are remarkable here: (1) A small perturbation influences very much the 
evolution of the system with t — > oo; in a sense it defines the character of 
asymtotic solutions. (2) One can say that in spite of the fact that a small 
perturbation is a rapidly damping one, the stable asymptotic solution 
"remembers" the initial perturbation. In other words, the nonlinear 
system, i.e., a LJJ, described with the sine-Gordon equation shows an 
effect of memory. However, in Ref. 12, the LJJ is studied solely under 
the influence of an external constant magnetic field. Therefore, it is of 
interest to investigate the LJJ from the point of view of the effect of 
memory not only in the presence of an external constant magnetic field 
but also under the influence of a dc bias current through the junction 
causing an exitation of dynamical chaos. How will the effect of memory 
in the presence of a dc bias current be shown? Will this effect take place 
in the states of dynamical chaos in general? Below we will try to give 
answers on these questions. 

We write down the sine-Gordon equation in the presence of a dc bias 
current in a LJJ in the form 



where ip(x, t) is the Josephson phase variable, x is the distance along the 
junction normalized to the Josephson penetration length Aj, 



$o is the flux quantum, j c is the critical current density of the Josephson 
junction, d = 2A^ + 6, A^ is the London penetration length, b is the 
thickness of the dielectric barrier, t is the time normalized to the inverse 
of the Josephson plasma frequency uj, 



C is the junction capacitance per unit area, 7 is the dissipative coefficient 
per unit area, and (3 is the dc bias current density normalized to j c . 



tp tt (x, t) + 2j(f t (x, t) - (p xx (x, t) = - sin ip(x, t) + f3, 



(1) 
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We write down the boundary condition for Eq. (1) in the form 

d -^\ x= o = H(0,t) = d -t^\ x =L = H(L,t) () 
= Ho(l-ae-^ to coat), { ' 

where L is the total length of the junction normalized to Xj, H is 
the external constant magnetic field perpendicular to the junction and 
normalized as well as H(0,t) and H(L,t) to 2 ^\jd ■> a is the controlling 
(perturbation) parameter characterizing the rapidly damping in the time 
perturbation, and to is the characteristic time of this damping perturba- 
tion normalized to uj 1 . 

Eq. (1) with boundary condition (2) is solved numerically. In con- 
trast to the case (3 = considered by us 12 , the picture of magnetic field 
evolution in the junction turns out more complicated by (3 ^ as will be 
shown below. Physically this is connected with the fact that the energy 
balance in the junction in the presence of bias current is such that the 
energy brought into the system with this current can make up for the en- 
ergy loss because of dissipation or exceed it. And so one can expect that 
at few values of (3 a state of junction will differ little from the station- 
ary one described by the Ferrell-Prange equation. At sufficiently large 
values of (3 asymptotic regular (periodic) solutions that represent waves 
- fluxons and antifluxons — moving along the junction and interacting 
among themselves and with junction boundaries and also nonregular so- 
lutions that represent dynamical chaos will take place 1 . Our calculations 
showed that if an asymtotic state at a = (further we shall call the state 
at a = as "starting") is regular, by introduction at the initial moment 
of a rapidly damping perturbation defined with the controlling parame- 
ter a as well as in a stationary case, examined in Ref. 12, the selection 
of the asymtotic solution by a given set of parameters Hq, 7, L, and (3 
is determined with this parameter a; i.e., the system "remembers" the 
form of the rapidly damping perturbation and chooses the way of further 
evolution in accordance with this. (It is noteworthy that the effect of 
memory discussed here happens in a dissipative system and so it is not 
connected with the reproduction of a signal as it takes place, for exam- 
ple, in noncollision plasma in the effect of plasma echo 14 ). Furthermore, 
our calculations show that if the "starting" state is regular (periodic), 
stationary states can also arise by introduction of perturbation (a 7^ 0) 
which it is of surprise in itself. However, the most remarkable fact is 



3 



that the system is very sensitive to a rapidly damping perturbation as 
well as in the dynamical chaos conditions; i.e., the system has memory 
in this case too. 

For the quantitative description of different characteristic states we 
used Lyapunov indices. We write down Eq. (1) in the form 



<Pt = v, 

V t = -2-/V + ip xx - sin ip + ft, 



(3) 



or, that is the same, in the form 

* = F(z), (4) 

where z is the vector with the components if and V, z = f f J , and 



V 



F(z) is defined as follows: 

F(z) = F(tp,V) = 



V 

-27 + (f xx - sinif + (3 



Let z{t) be a solution of Eq. (4). Then we can write down the equation 
for variations: 

dF(z(t)) 

w t = -^w, 



w — 



wi \ (5) 

W 2 



We define the Lyapunov index (LI) as 

A=limilnM, (6) 

t-°o* IKo)||' v 1 

where ||w;|| is the vector norm that we define as Euclidean norm 

L 

\\w\\ 2 = Uw\ + w\)dx. (7) 



o 



Depending on direction of the initial vector w(0), different Li's will exist 
and their number will be infinite. The definition of LI (6) for system 
(3) is a natural generalization of a LI for finite-dimensional dynamic 
systems. 13 

The maximum LI plays a very important role because it is precisely 
this maximum that determines the motion character - exponential 



growth, decay, or zero change - - for the majority of the trajectories 
of the system. A set of initial data w(0) for which formula (6) gives Li's 
differing from the maximum one is negligibly small and by numerical 
calculations this formula gives the maximum LI as a rule. 

The results a numerical solution of Eq. (1) with boundary conditions 
(2) and calculation of LI (6) showed that there exist three forms of 
characteristic states of the system for which the maximum LI can be 
as follows: (1)A > 0, (2)A < 0, and (3)A < 0. The states with A > 
represent the dynamical chaos states (Fig. 1), the states with A < 
represent the stable stationary states (Fig. 2), and the states with A < 
represent the regular (periodic) states (Fig. 3). All states in Figs. 1-3 
that are shown as illustration are as "starting" ones (a = 0) and they 
are calculated with identical values of the parameters Hq = 1.25, L = 5, 
and 7 = 0.26, but with different values of (3. For the chaos state in Fig. 
1, j3 = 0.50; for the stationary one in Fig. 2, j3 = 0.427; and for the 
regular one in Fig. 3, (3 = 0.60. In Figs. 1-3 the dependences of the 
potential (ft (potentials are normalized to the value V c = ^) on time 
and the calculated values LI A corresponding to them are shown. We 
note that the "starting" chaotic state, represented in Fig. 1, is the same 
as in Ref. 1. 

Let us examine the specific set of parameters, i.e., the definite point 
in parameter space, corresponding to the chaotic "starting" state: Hq = 
1.25, L = 5, 7 = 0.26, (3 = 0.44, and a = 0. If we input now the 
rapidly damping in the time perturbation determined by the control- 
ling parameter a, the system does not remain in the previous chaotic 
state as the calculations show, but wanders between all three forms of 
states: chaotic, stationary, and regular, when this parameter changes. 
By calculations the following hierarchy of times was used: to T r T, 
where r r is the characteristic time of relaxation processes (r r is the time 
of relaxation to asymtotic states) and T is the time of observation. The 
values of characteristic times in our calculations were as follows: to = 5, 
7> ~ 60, T = 2000. At first sight, one might have expected that initial 
perturbations at the time interval T damping at the time about to would 
be forgotten and they will have no influence on the evolution in the large 
time interval. [We notice that r r is not equal to 7 _1 , exactly r r ^> 7 _1 in 
our case. This is connected with the system's nonlinearity. The value of 
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T r is defined from numerical calculation of our problem (1), (2).] In Fig. 
4 are shown the results of calculation of Li's with values of parameters 
for Hq, L, 7, and (3 mentioned above and by specific values of param- 
eter a. As we can see from Fig. 4 three typical clusters of states take 
place: a cluster of chaotic states ch (in Fig. 4 the following values of the 
parameter a correspond to them: a =0, 0.175, 0.180, 0.280), a cluster 
of regular states r (a = 0.290, 0.300, 0.320), and a cluster of stationary 
states s (a = 0.190, 0.195, 0.285). For the cluster ch the values of the 
LI A « 5 • 10~ 2 ; for the cluster r, A « — 10~ 3 ; and for the cluster s, 
A « -lO -1 . 

In Fig. 5 the potentials on the junction (p t depending on time for three 
values of parameter a differing from each other on 0.005 and belonging 
to cluster ch, a = 0.280 [Fig. 5a]; to cluster s, a = 0.285 [Fig. 5b]; and 
to cluster r, a = 0.290 (Fig. 5c) are showed; their Li's are represented 
in Fig. 4. Thus, a small change of parameter a leads to a transition 
between all three characteristic states of the system. 

In Table I transitions between chaotic ch, stationary s, and regular 
r states are represented by a change of the parameter a from 4.000 to 
4.155 at the fixed remaining parameters indicated above. 

Table I. States of a LJJ 
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We note that the transitions between the states, reduced in Table I 
and stipulated by a change of the perturbation parameter a, correspond 
to the "starting" state of chaos (a = 0). Thus, by the given values of 
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Hq, L, 7, and (3 a final asymptotic state is determined with the pa- 
rameter a independently of this, whether the "starting" state is chaotic 
or not. Such asymptotic behavior of the system says that dynamical 
chaos essentially differed from statistical chaos by which any perturba- 
tion damps rapidly and a system relaxes to its final state (for example, 
to the state of thermodynamical equilibrium) completely "having forgot- 
ten" an initial perturbation; i.e., the final state does not depend on this 
perturbation. As we see, a system in the state of dynamical chaos in 
contrast to a system being in the state of statistical chaos "remembers" 
the initial perturbation and, in any sense, final states and transitions 
between them are defined with the very initial perturbation. It makes it 
possible for us to recognize this memory effect in the system described 
with the sine-Gordon equation with dissipation and in the presence of an 
external magnetic field and a bias current. Thus, the dynamical chaos 
originating in a nonlinear system does not destroy an initial information; 
i.e., the nonlinear system has a memory in the states of dynamical chaos 
as well. 
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1. The potential <p t (a) and the Lyapunov index (b) in a chaotic state 
by (3 = 0.50. The values of other parameters are Hq = 1.25, L = 5, 
7 = 0.26, a = 0. 

2. The potential ip t (a) and the Lyapunov index (b) in a stationary 
state by (3 = 0.427. Other parameters are the same as those in Fig. 
1. 

3. The potential ip t (a) and the Lyapunov index (b) in a regular state 
by (3 = 0.60. Other parameters are the same as those in Fig. 1. 

4. The Lyapunov indices A vs parameter a: ch is the cluster of chaotic 
states (a =0, 0.175, 0.180, 0.280), r is the cluster of regular states 
(a =0.290, 0.300, 0.320), and s is the cluster of stationary states 
(a =0.190, 0.195, 0.285). The values of other parameters are just 
the same: H = 1.25, L = 5, 7 = 0.26, (3 = 0.44. 

5. The potential <p t vs t belonging to the cluster ch, a = 0.280 (a), to 
the cluster 5, a = 0.285 (b), and to the cluster r, a = 0.290 (c). 
Other parameters are the same as those in Fig. 4. 
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